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1 Introduction 



The Standard Model (SM) of particle physics is very successful in explaining most of the elemen- 
tary particle physics phenomena at the electroweak scale. Unfortunately, despite of its eminent 
success, the SM suffers from certain drawbacks, the hierarchy problem being a major one. Su- 
persymmetry (SUSY) ^ provides an elegant solution and, consequently has been extensively 
studied as a model beyond the SM. If such a new physics exists at the electroweak scale, then 
it may provide some experimental signature in the future, or even in present data. Looking for 
such new physics effects constitutes a major area of research in high energy physics today. There 
are two major ways one can see such effects. One is to observe their direct effects in the high 
energy collider experiments when a new particle is produced and observed through its decays. 
The other way, is to look for indirect evidence in the deviation from the SM prediction of low 
energy experimental data. 

Recent results from the ongoing experiments in i?-physics at CLEO, BaBar, and Belle have 
attracted lot of attention. Much of this attention has been devoted to the results on nonleptonic 
decays of B mesons, which can be used to extract information on the CKM matrix elements 
and CP violation. The theoretical understanding of the nonleptonic decays of B mesons is an 
extremely demanding challenge due to difficulties in calculating the relevant hadronic matrix ele- 
ments. To have some idea of the magnitude of the matrix elements, one usually uses factorization 
method, factorizing the four quark operators relevant to non-leptonic B decays into the products 
of two currents and evaluating separately the matrix elements of the two currents. Recently the 
QCD improved factorization method for the hadronic B decays has been developed. This method 
incorporates elements of the naive factorization approach (as its leading term) and perturbative 
QCD corrections (as subleading contributions) allowing one to compute systematic radiative cor- 
rections to the naive factorization for the hadronic B decays ^ . In our analysis we will use the 
formalism developed in Ref. Q. This QCD-improved factorization method improves the analysis 
on several aspects, including among others the number of colors, the gluon virtuality, the renor- 
malization scale, and the scheme dependence. The method is expected to give a good estimate of 
the magnitudes of the hadronic matrix elements in non-leptonic B decays, and has been used to 
calculate B decays in the SM [Q, ^, ^, ^ and models beyond ||7j. 

The construction of most general supersymmetric extension of the standard model leads to 
baryon (B) and lepton (L) number violating operators in the superpotential. The simultaneous 
presence of both (L) and (B) number violating operators induce rapid proton decay which may 
contradict the strict experimental bound In order to keep the proton lifetime within the 
experimental limit, one needs to impose additional symmetry beyond the SM gauge symmetry to 
force the unwanted baryon and lepton number violating interactions to vanish. In most cases, this 
has been done by imposing a discrete symmetry called i?-parity defined as i? = (— i)3B+l+2S^ 



where, S is the spin of the particle. This symmetry not only forbids rapid proton decay ||l^, it 
also render stable the lightest supersymmetric particle (LSP). However, this symmetry is ad hoc 
in nature. There are no strong theoretical arguments in support of this discrete symmetry. Hence, 
it is interesting to see the phenomenological consequences of the breaking of i?-parity in such a 
way that either B or L number is violated both are not simultaneously violated, thus avoiding 
rapid proton decays. Extensive studies have been done to look for the direct as well as indirect 
evidence of i?-parity violation from different processes and to put constraints on various i?-parity 



violating couphngs |lT|, |12[|2C|. 

The main purpose of this paper is to constraint various i?-parity violating couplings using the 
data on i? — > PP decay channels and a calculation based on QCD-improved factorization. Here 
P is one of the S(3) flavor octet pseudoscalars. We find that using the experimental data on the 
branching ratios of ^ PP mode, stringent upper bounds on the products of several L and B 
violating couplings can be obtained. Many of the bounds obtained are stronger than the existing 
ones. 

The organization of the paper is the following. In section (2) we study possible four quark 
operators which can induce B PP decays with i?-parity violating interactions. In section (3) 
we calculate B — > vrvr. Kit, KK decay amplitudes using the QCD-improved factorization method. 
At last in section (4), we carry out numerical analysis, present and discuss our results and finally 
draw our conclusions. 



2 New operators for B PP with i?-parity violation 

The most general superpotential of the Minimal Supersymmetric Standard Model (MSSM) Q, 
which describes SU (3) x SU{2) x U{1) gauge invariant, renormalizable and supersymmetric theory, 
with minimal particle content, has -R-parity violating interaction termsjP]: 




where, L and Q are the 5'[/(2)-doublet lepton and quark superfields and E'^, and D'^ are the 
singlet superfields, while A; are flavor indices. In writing the above we have omitted gauge 
indices, which ensures that the Xijk are antisymmetric in i and j, and the A" are antisymmetric 
in j and k. It is very clear that the first two terms in Equation (||) violate lepton number, while 
the last term violates baryon number. 

The above -R-parity violating interaction, in general can have 27 A'-type and 9 each of A and 
A"-type of new couplings. Not all of them will induce B — > PP decays to the lowest order. At the 
dimension six four-fermion interaction level only operators corresponding to the couplings A- 



and X'l^j^ will lead to hadronic B PP decays. They are given by 
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where is the sfermion mass, L{R) = (1 T Ib)/'^^ and a, (3 are the color indices. 
There are three types of four-quark operator in the above four quark interactions, 

Oa = iPaq/3)v±A{r^ba)vTA- 
Ob = {paqa)v+A{ri3bi3)v+A, 

Oc = {Paqi3)v+A{fab/3)v+A, (3) 

where {pq)v±A = Plix{^ ± Ts)^- 

The first operator above occurs in A' interactions and the last two are in A" interactions. The 
above operators are evaluated at the common sfermion mass scale {mj) of 100 GeV. At a scale 
around /x = mi,, these operators will induce nonzero matrix elements causing B — ^ PP decays. We 
denote c{^)a,b,c the Wilson coefficients of the operators Oa,b,c at the scale /x. Renormalization 
group running of these coefficients from to /x will modify them. For c{ij.)a, we have 

c(m)a = r?-8/*c(m.)A, (4) 



where 77 = as{mj)/as{n), and /5o = 11 — 2//3 with / the number of quarks with mass below fi. 
The other two coefficients c{mj)B and c{mj)c when evolved down to the scale from the high 
scale m j will mix and are given by 

c(m)b = l[v^^^"{c{mj)B + c(m^-)c) + r?-^/^°(c(m^-)B - c(m^-)c)], 

cifih = \[v''^^"ic{mj)B + c{mj)c) - r^/^'\c{m^)B - c{m^)c)]. (5) 

To obtain the B PP decay amplitudes induced by these operators, one needs to evaluate the 
related hadronic matrix elements. In the next section, we will use the QCD-improved factorization 
method to carry out the analysis. 



3 B ^ PP decay amplitude with i?-parity violation 

The factorization approximation has been used to provide estimates for the decay amplitudes in 
B decays. Recently it has been shown that factorization approximation in fact is supported by 
perturbative QCD calculations in the heavy quark limit, and the QCD-improved factorization 
method has been developed [||, Q|. This new factorization formula incorporates elements of the 
naive factorization approach and introduces corrections in the amplitudes. 

In the heavy quark limit, the decay amplitude B ^ P1P2 due to some particular operator Oi 
can be represented in the form 



< P1P2 \ Oi\B >=< P2 I J2 I >< Pi I Ji I 5 > 



1 + ^ r„< + ©(AqcdM) 



(6) 



The above result reduces to the naive factorization if we neglect the power corrections in Agc/j/mb 
and the radiative corrections in Ug- The radiative corrections, which are dominated by the hard 
gluon exchange can be computed with perturbation theory in the heavy quark limit, in terms of 
the convolution of the hard scattering kernel and the light cone distribution amplitudes of the 
mesons. Then a factorization formula for B P1P2 decay can be written as 

<PiP2\Oi\B> = F^-^^'iO) f\xTl{x)<^p,{x) (7) 
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In the above formula, ^BiO and ^p.{x) (i = 1,2) are the leading twist light cone distribution 



amplitudes of the S-meson and the light pseudoscalar mesons |21, ^] respectively, and the T- ' 
describes the hard scattering kernel which can be calculated in perturbative QCD g |, |]. The 
diagrams generating the hard scattering kernels T^'^^ in the SM are shown in Figure 1. Figures 
l(a)-(d) depicts vertex corrections. Figures 1(e) and 1(f) penguin corrections, and Figures 1(g) 
and 1(h) hard spectator scattering. 



Applying the QCD-improved factorization method to the operators in the previous section, 
we obtain the contributions from i?-parity violating interactions. Since we are considering the 
leading effects, we need only evaluate Figures l(a)-(d) for the vertex corrections and Figures 1(g) 
and 1(h) for the hard-spectator scatterings. The penguin types are higher order corrections. The 
results are listed below: 
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where /j, Fq are decay constant and form factors, respectively. The parameters Ri and rj are 
given by, 
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with q = u for charged kaon and q = d for neutral kaon. The parameters a'- and a' are defined as 
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Here A'c = 3 is the number of colors, Cp = {N^ - 1) /2Nc, Pi is the final state meson absorbing the 

light spectator quark from i?-meson, while P2 is another final state light meson which is composed 

of the quarks produced from the weak decay of b quark. 
(') 

Vp comes from vertex corrections {P = it, X) (first four diagrams of Figure 1) is give by 



Vp 
Vp 



12 In 



12 In 



3(- 



- 2x 



18 H 
6 + 



dxg{x)(j)p{x), 



dxg{l - x)(f>p{x), 



Inx — iir) + 



X 



2Lio(x) 



,9 21nx , , , 

In^ X H (3 + 2z7r) - (x ^ 1 

1 — X 



(23) 
(24) 
(25) 



The contributions from the hard spectator scattering as shown in the last two diagrams in 
Figure 2 give leading twist and chirally-enhanced twist-3 contributions to T/^ parametrized by 
hI^Ip^. The detailed expressions can be found in H. 

From the above, the amplitude Ap{B^ — > K^K^) is zero without annihilation contributions 
(shown in Figure 2). When annihilation contributions are included, this amplitude becomes 
nonzero and additional contributions to other decay amplitudes are generated. We list these 
contributions in the following: 
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The parameters Ai are given by 
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where R.,^, Xa = dy/y parameterizes the divergent endpoint integrals, we take the value 
Xa = ln(mB/A/j) with = 0.5 GeV being the typical hadronic scale. 



4 Results and Discussion 

In this section we discuss the constraints on the product of a pair of i?-parity violating couplings. 
With i?-parity violating interactions, the total B PP decay amplitudes contain the SM part 
of the decay amplitudes Asm which have been obtained in Ref. Q plus the R-parity violating 
amplitudes obtained in the previous section. Comparing the branching ratios obtained theoret- 
ically with known experimental data, constraints on the i?-parity violating interactions can be 
obtained. 



As we have shown in the previous section that several combinations of -R-parity violating 
couplings can contribute to a particular charmless hadronic B PP decay. In our numerical 
analysis, we assume that only one pair of ii-parity violating couplings are nonzero at a time. 
This restriction may seem to be unnatural, however, it is an useful approach that allows one a 
quantitative feeling of the various experimental constraints. 

For the numerical computations we use the averaged value of the experimental input as shown 
m Table |l[ We average the data from BaBar, Belle, CLEO B factories as their results are 
uncorrelated. If the experimental branching ratio has only upper limits, we select the most 
stringent one. It is clear from the Table |l] that for B vt^tt'^, B tt^tt^ and B KK modes 
have only the 90% C.L. upper limit while others have results at 95% C.L. 

In our numerical calculation we use mf,(mfe) = 4.2 GeV for the b quark mass. For the 
other lighter quark masses we take their central values mdmi,) = 1.3 it 0.2 GeV, ms{2 GeV) = 
110.0 ± 0.25 MeV, (m„ + md)(2GeV) = 9.1 ± 2.1 MeV. For the decay constants g and form 
factors||, we use = 0.131 GeV, fx = 0.160 GeV, fs = 0.180 GeV, F^^"" = 0.28, 
^ttK — pB^wj^ = 0.9. For the KM matrix elements, we fix A = 0.2196, and take the central 
values of | Kb |= 0.0402 ± 0.0019 | ^ |= 0.090 ±0.025 g, but allow the CP violating phase 7 to 



vary from to 2tt. In the SM, the CP violating phase 7 is well constrained However, with 
i?-parity violation this constraint may be relaxed. To accommodate this we vary 7 from to 2tt 
to obtain conservative limits on i?-parity violating couplings. 

We also need to know various light cone distribution amplitudes. The leading-twist light cone 
distribution amplitudes of the light pseudoscalar mesons [^] can be expanded in Gegenbauer 
polynomials. We truncate this expansion at n = 2. 

(t>p{x,fi) =6x{l-x) [l + af(/i)C7f/^^(2x-l)+a|'(/u)cf/^^(2x-l)] , (35) 

where cj'^^^^(u) = 3u and C^^'^\u) = |(5n^ — 1). The distribution amplitude parameters 

for P = K,TT are: of = 0.3, af = 0.1, aj = and = 0.1 §. 

The hard spectator contributions to the coefficients a" 2 , and a' are parametrized in terms of 

(') 

a single (complex) quantity Hp^p^, which suffers from large theoretical uncertainties related to 
the regularization of the divergent endpoint integral. Following [0] we use Ht^k = Hkk = 0.99 
at the scale /i = rrih and 

H'p^Pi = -ffp2Pn -f^TTTT = Hk-k = rnKH-irK- (36) 

Here we should like to discuss the numerical changes of the coefficients a'( 2 , and a' due to our 
use of QCD improved factorization calculations instead of the naive factorization approximation. 
For example, in the B vrvr decay mode using the naive factorization scheme the coefficients a'/, 
and a' have values 0.939455 and 0.661932 respectively. Using the QCD improved factorization 
method we get 1.0164 + 0.105414 i and 0.581234 + 0.148545i for a'{ and a' respectively. For the 



coefficient a'/ we get an enhancement of the order 10% in the real part and a new imaginary 
contribution, while for the coefficient a' we have a negative contribution of the order 12% in 
the real part but a new imaginary contribution as before. The QCD improved factorization 
gets contributions from vertex corrections, hard scattering and from weak annihilation. The 
contribution from annihilation diagrams is small compared with the current-current contributions. 
It is less than 5% in decays involving a'/ 2 while for decays only involving a' the annihilation 
contributions can be as large as 20%. 

The results are shown in Tables ^ and ^ In the Table ^ we display limits on different pair 
of A"-type couplings at 95% C.L. For the published experimental data with 90%C.L. it is hard to 
get the upper limit corresponding to 95%C.L., so we denote those bounds by (*) in the Table |3|. 
In Table |^ we show the bounds for the product of A'-type of couplings. The bounds are obtained 
assuming 100 GeV common sfermion masses, so for other values of the sfermino mass, the bounds 
on the couplings can be easily obtained by scaling them by m^/(100GeV)^. 

We note that for several decays i?-parity violation contributes only through the annihilation 
contribution. There are possible large uncertainties in the evaluation of the weak annihilation 
contribution to charmless B decays due to the fact that they are power suppressed in the heavy 
quark limit and the weak annihilation effects exhibit endpoint singularities. We therefore denote 
the modes which only receive annihilation contributions with (f) to remind that there may be 
large uncertainties there. 

From the Table ^ and Table Q we choose the pairs of i?-parity violating couplings with most 
stringent bounds and we display them in Table |5| with the existing limits on such product of 
couplings. Among the existing limits, there are several cases, where there is no direct limit on 
the products of the couplings, in that case we take the products of their individual bounds from 
Ref . ||Tl|| . As it can be seen in [11| that the bounds on the first two generation individual couplings 
are much stronger than the third generation. For this reason, in most of the cases, the product 
of either or both first two generation individual couplings are much stronger than our prediction. 
Furthermore, we find that bounds on the pair of couplings A'/^2'^ii3) ^"^^ ^112-^123 obtained from 
n — n oscillation and double nucleon decays, A^3^A-22, '^232^^211 ^332-^3115 obtained from Anix 
are stronger than ours. It is, however, interesting to note that, more than half of the bounds 
obtained in this paper are better than the existing ones. In addition to that, the bounds on 
-^212^2131 ^312-^313' -^212^223' ^312 ^^323 Couplings are completely new, in a sense that there 
were no previous bounds on these pairs of couplings from any experimental data, they were from 
perturbative unitarity. In two of these cases (A212A213 and X'^i2X'^i^) we have improved the existing 
bound by over two orders of magnitude. 

Before we go conclude we would like mention the possible theoretical uncertainties in the 
bounds arising from several input parameters. As described before, we have obtained all of our 
bounds by using the central values for the relevant parameters, the CKM matrix elements, the 



quark masses, the form factors. We have also parametrized the hard spectator contribution by 
some default number. Some of these parameters may change up to 20%, and the uncertainty in 
the hard spectator contribution may be even larger. Even allowing for this, we do not expect our 
result to change more than factor of two. 

In this paper we have studied B ^ tttt, B ^ ttK and B KK decays based on the QCD 
improved factorization approach in the presence of i?-parity violating couplings. Comparing our 
calculated branching ratios B PP with the experimental data, we have obtained bounds on 
the products of two i?-parity violating couplings. We found that most of the bounds we obtained 
on combinations of the A' and A" type of couplings are stronger than the existing limits. Rare 
hadronic B decays can provide important information about i2-parity violating interactions. 
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Table 1: Experimental results for the CP-averaged B 
ratios in units of 10~^|25|. 



vrvr, B ttK and B KK branching 





B TTTT 


B ^ Kit 


B^ttK 


B KK 


c{mb)A 


1.986 


c{mb)B 


1.409 


c(mb)c 


-1.409 


a'l 


1.016 + 0.105i 


1.002 + 0.105i 


0.998 + 0.137i 


1.020 + 0.074't 


n 
O2 


-1.016 -0.105i 


-1.002 -O.lOSi 


-0.998 - 0.137i 


-1.020 - 0.074i 


a' 


0.581 + 0.149« 


0.561 + 0.149i 


0.607 + 0.104i 


0.535 + 0.1931 




O.OiG 




-0.016 


rAb'3 


0.141 


TAb'i 


0.022 



Table 2: Wilson coefficients, and related parameters due to i?-parity violating interactions in 
different decays, /i = 7715, assuming 100 GeV sfermion mass. The first final state meson is made 
of spectator quark, while the second final state meson is formed of quark pairs originating from h 
quark weak decay vertex. In the above ta = fBf-jT/{iTT-B-^^~*^i^))- 
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Table 3: 95% C.L. limits on the products of A"-type i2-parity couplings assuming 100 GeV common 
sfermion mass. The limits correspond to the label (*) are at 90% C.L., while the limits denoted 
by (t) are obtained from decay modes with only annihilation contributions. 
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Table 4: 95% C.L. limits on the products of A' i?-parity violating couplings for a common sfermion 
mass mj = 100 GeV. The limits correspond to the label (*) are at 90% C.L., while the limits 
denoted denoted by (f) are from decay modes with only annihilation contributions. 





(JUi iillilLb 


jJUJVlULlb illlllLh 


process of others constraints 


'^112'^113 


[U.OU A ±U J 


9 V 1 C\~^ 

Zr A ±U 


n — n and double neucleon decay ]l9[ 


x" x" 


1 0'^ y 1 0""^ 

-L.UO A ±\J 


1 


perturbativity bound [ 


Li 


X X 

-^312^313 


1 ff^ V 1 0"^^ 

±.UO A ±\J 


1 ^ 


perturbativity bound [ 


i| 


x" x" 

'^112'^123 


[O.OU A lU J 


9 V 1 

Zr A ±U 


double nucleon decay | 


l| 


X X 

'X212-^223 


Q 1 Q V 1 n~^* 

o. ly A ±\j 


1 ^ 


perturbativity bound [ 




X X 

'^312^^323 


Q 1 Q y 1 f)~'-^* 

o. ly A ±u 


1 ^ 


perturbativity bound [ 




X X 

-^111-^113 


[i.'iy A lu J 


1 1 V 1 




\' x' 

'^211^213 


1 AQ V 1 n~^* 

A ±\J 


O.O A ±U 




X X 

'^311'^313 


1 /iQ V 1 n~^* 

L.^U A lU 


O.O X lU 


Am ^ IrXll 


x' x' 

-^111^123 


[Z.io A lU J 


9 9 V 1 




X X 

-^211-^223 


Zr. ±0 A ±\J 


1 9 V 1 




'^311'^323 


9 1 Q V 1 

Zr. ±0 A ±\J 


±.\J A ±U 


Am fxll 


X X 

-^lll'^131 


[i.yo A lu J 


1 V 1 

l.U A lU 




x' \' 

-^211'^231 


1 Qf^ y 1 0""^ 

I.yo A lu 


1 1 V 1 




X X 

'^311-^331 


1 Qf^ V 1 
i.yo A lu 


Pi y 1 
O.U X lU 




X X 
-^111^132 


[9 1 s V 1 n~^l 


1 /I \^ 1 0""^ 




X X 

-^211 ^^232 


[9 1 Q V 1 n~^l 

[Z.±0 A lU J 


ZL 7 y 1 0""^ 


Am I.' ll^rll 


x' \' 

-^311 ^^332 


[Z.io A lU J 


/I 7 \^ 1 0""^ 


A m ^ - [brii 


X X 

'^112-^113 


[1 . ( i A lU J 


/I /I y 1 n~*^ 




x' x' 

'X212'^213 


1 71 V 1 

1. / _L A lU 


^ \^ 1 

O.O X lU 




\ X 

'^312-^313 


1 71 1 


i.z X iU 




x' x' 

'^112'^132 


70 y 1 O"'^* 
D. / U A lu 


Pi Q y 1 

o.y X lu 




'^212^^232 


(S 70 V 1 0""^* 


Q q V 1 n~-^ 

O.O A ±U 




-^312'^332 




5.0 X 10^^ 




'^113'^122 


[9 y 1 0^^*1 


Q y 1 n~'^ 




'^213'^222 


9 y 1 n~^* 


1.2 X 10~^ 




-^313-^322 




5.7 X 10~^ 




-^121-^123 


6.70 X 10"''* 


1.4 X 10"^ 




■^221-^223 


6.70 X 10"''* 


1.4 X 10"^ 




-^321-^323 


6.70 X 10"''* 


1.4 X 10"^ 




■^121-^131 


[1.33 X 10"^] 


8.2 X 10"-^ 




-^221-^231 


1.33 X 10"=* 


3.2 X 10"^ 




■^321-^331 


1.33 X 10"^ 


2.3 X 10"^ 




-^122-^131 


[2.88 X 10"^*] 


1.0 X IQ-^ 


AmK@ 


■^222-^231 


[2.88 X 10"^*] 


1.0 X IQ-'* 


AmK||l 


-^322-^331 


[2.88 X 10"^*] 


1.0 X IQ-^ 


AmA-|c|] 



Table 5: Comparison of constraints obtained in this paper with other existing constraints. 95% 
C.L. limits on the products of different -R-parity couplings for sfermion masses "ij = 100 GeV. 
The limits correspond to the label (*) are at 90% C.L.. The limits shown in the square bracket 
are weaker than the existing bounds. The bounds corresponding to the processes not shown in 
the Table are coming from the product of individual couplings [p!9|]. 
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Figure 1: Order of Os corrections to the hard- scattering kernels and T^^ . The quark lines 
directed upwards represent the ejected quark pairs from weak decays of b-quark. 
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Figure 2: Annihilation diagrams 



